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ABSTRACT 

The usual quantizer based on an n-dimensional lattice A maps a point x G R n to a closest 
C — ' 

| lattice point. Suppose A is the intersection of lattices Ai, . . . , A r . Then one may instead combine 

. the information obtained by simultaneously quantizing x with respect to each of the Aj. This 

i -pH | corresponds to decomposing W 1 into a honeycomb of cells which are the intersections of the Voronoi 

cells for the Aj, and identifying the cell to which x belongs. This paper shows how to write several 
standard lattices (the face-centered and body-centered cubic lattices, the root lattices D4, Eq, E$, 

•i-H . 

^ ■ the Coxeter-Todd, Barnes-Wall and Leech lattices, etc.) in a canonical way as intersections of a 

b : 

small number of simpler, decomposable, lattices. The cells of the honeycombs are given explicitly 
and the mean squared quantizing error calculated in the cases when the intersection lattice is the 
face-centered or body-centered cubic lattice or the lattice D4. 

1. Introduction 

An n-dimensional lattice A determines a partition of M. n into cells which are congruent copies 



of the Voronoi cell containing the origin. A vector quantizer based on A ( plfl , [16], [|18|]) is a map 



Qa : M. n — > W 1 which maps x £ W 1 to the closest lattice point, i.e. to the lattice point at the center 



"This work was carried out during B. Beferull-Lozano's visit to AT&T Labs in the summer of 2001. 
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of the Voronoi cell containing x. (In the case of a tie, one of the closest lattice points is chosen at 
random.) 

Suppose A is the intersection of n-dimensional lattices Ai, . . . , A r . A recent paper j3j discusses 
(among other things) the possibility of replacing the quantizer Q\ by a "multiple description 
quantizer" which simultaneously quantizes with respect to each of the Aj, i.e. computes 



(Q\i (x) , Qa 2 (x) , . . . , Q\ r (x) ) . 



(1) 



This gives rise to a different partition of M. n : the cells are now the intersections of the Voronoi 
cells of the individual Aj, and (|l]) specifies the cell to which x belongs. 
There are several reasons for investigating quantizers of this type. 

(i) If the A, are simpler lattices than A then (|l]) may be easier to compute than Q\(x). 

(ii) Perhaps this approach will lead to new insights on the so-far intractable problem of finding 
good lattice quantizers in high dimensions (cf. [jTT] , Chap. 2]). Even in 24 dimensions the best lattice 
quantizer presently known, the Leech lattice, is quite complicated to analyze and to implement - 
its Voronoi cell has 16969680 faces and over 10 21 vertices (0, Chaps. 21, 22, 23, 25], §, @, ||). 



(iii) The individual Q\ i (x) could be communicated over separate channels; in the event of one 
or more channels failing a reasonably good approximation to x will still be obtained. Other multiple 
description quantizers have recently been studied in, for example, [O], |2S]. 



(iv) Initially it seemed possible that this approach might lead to quantizers with a lower mean 
squared error that can be obtained from lattice quantizers. The investigations reported here now 
make this unlikely (but see the remarks in the final section). 

However, our focus here is not on applications but on the theoretical aspects of these quantizers, 
and in particular on two questions: which lattices have a "nice" description as the intersection of 
simpler lattices, and what are the associated "honeycombs," or decompositions of space into cells? 

Reference J|] gives one very appealing example, shown here in Figure [l]. The familiar planar 
hexagonal lattice A2 (large circles) can be obtained as the intersection of three rectangular lattices, 
all rotations of each other by 120°: these are the lattices generated respectively by the two vectors 
OA, the two vectors OB and the two vectors OC. If we use generator matrices to specify these 
lattices (the rows span the lattices) then the three rectangular lattices Ai, A2, A3 have generator 
matrices 

r./s ni VE l 1 r ^3 1 " 

(2) 
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and their intersection has generator matrix 



which is indeed a copy of the A2 lattice. 
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Figure 1: The hexagonal lattice (heavy circles) as the intersection of three rectangular lattices 
(spanned by the vectors OA, OB and OC resp.). 

In this case the honeycomb contains four kinds of cells, as shown by the heavy (solid) lines in 
Fig. ||. For example, if the point x being quantized is close to (0,0), then Q\ 1 (x) = Qa 2 ( x ) = 
Qa 3 (x) = (0, 0), and the cell containing x is the intersection of the Voronoi cells containing (0, 0) of 
the three Aj. This is the horizontally shaded hexagon in Fig. ||. Just to the North of this hexagon 
the cell is the intersection of the Voronoi cell for Ai that contains (0, 1) with the Voronoi cells at 
(0,0) for A2 and A3: this is the small cross-hatched equilateral triangle. There are two further 
cells that are obtained in a similar manner: the diagonally shaded isosceles triangle and the larger 
vertically shaded equilateral triangle. This example will be discussed further in Section 4. 

The present paper was prompted by the question: how can Figures [I] and ^ be generalized to 
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Figure 2: Honeycomb associated with Voronoi cells in Fig. H. [Note to copy editor: Figures Hand 
H should be scaled so that the distance from to the point labeled A in Fig. ^ is the same as that 
distance in Fig. |].] 



higher dimensions? 

Section 2 gives some terminology and defines the normalized mean squared error G that will 
be used when comparing quantizers. Section 3 describes some general constructions for writing a 
lattice as an intersection and gives a number of examples, including the body-centered cubic (bcc) 
and face-centered cubic (fee) lattices and D3, the root lattices -D4, E§, E%, the Coxeter-Todd 
lattice K12, the Barnes- Wall lattices BW n and the Leech lattice A24. We focused attention on these 
lattices because A%, D* i: D4, Eq, Eg, K12, BW\q and A24 are the best quantizers currently known[] 
in their dimensions [11]. In fact A2 is optimal among all two-dimensional quantizers [O], and D% 
is optimal among three-dimensional lattice quantizers [||. 

In Section 4 we reexamine the quantizer of Figures [l] and |2[ Then Sections 5, 6 and 7 determine 
the honeycombs and mean squared errors for the bcc, fee and D4 lattices, respectively. The final 
section contains some conclusions and comments. 

Table [l] summarizes the main decompositions mentioned in this paper. 

We repeat that we do not propose that the quantizers described in this paper will have practical 
applications, nor that they are simpler than the usual lattice quantizers based on the intersection 
lattices. Our interest is in the geometry of the honeycombs associated with the quantizers. These 
do not seem to have been studied before^] — they are not mentioned for example in Okabe, Boots 

^Assuming always that the variable to be quantized is uniformly distributed over a large ball in R n . 
*Of course there is an extensive literature dealing with the Voronoi and Delaunay cell decompositions associated 
with various lattices — see for example Q , || , juj , Jnj , . 
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Lattice 


Copies 


Component 
lattice 


sections 


A 2 


3 


rectangular 


1,3,4 


A% (bcc) 


3 


"rectangular" 


3,5 


A 3 (fee) 


4 


"prismatic" 


3,6 


D 4 


3 


Z 4 


3,7 


E* 


4 


4* 


3 


Eg 


15 


Z 


3 


Es 


10 


A A 

A 2 


3 


Es 


5 


D\ 


3 


K12 


21 


A% 


3 


Leech 


4095 


Z 24 


3 


BW n 


nr=i 1 (2 j - 1) 


z n 


3 



Table 1: Summary of decompositions described in this paper. 



and Sugihara 23], nor in Wells [33], even in the chapter on Interpenetrating three-dimensional 
nets. There has been considerable interest among physicists in recent years in soap froths, both 
in connection with the study of colloids, dendritic polymers, etc., and in the Kelvin problem of 
finding the minimal-area honeycomb in R 3 ( pql , [p2[ , [pq] ). So it may be of interest to see a class 
of honeycombs that have a purely geometrical construction. 

We end this section by describing how the cells of the honeycombs in Sections 5-7 were found. 
Many cells could be found by elementary geometrical reasoning. But in complicated cases we 
carried out some or all of the following steps. 

(i) To find the cell containing a point x G R 3 or R 4 we first quantized x using each of the lattices 
Aj in turn. For each Aj, we determined the Voronoi cell Ci containing x, or, more precisely, the 
equations to the hyperplanes bounding Cj. 

(ii) Linear programming (in MATLAB) was then used to determine the vertices of the cell 
containing x. We let w range over a set of 130 points on a sphere centered at x (taken from the 
tables of spherical codes in |l^]), and, for each w, we maximized the inner product (w,z) subject 
to the constraints that z lie in the polytope formed by the intersection of all the hyperplanes found 
in (i). Any such solution z is a vertex of the cell, and since the w's are essentially random, the 130 
solutions should include all the vertices of the cell. 

(iii) The convex hull program QHULL (T^j was used to find the convex hull of these vertices. 
At this point we have candidates for all the cells in the honeycomb. Since it is theoretically 

possible (although unlikely) that step (i) might have failed to find all the vertices of a cell, we now 
verified by hand that the cells fitted together to form a proper tiling of the space. 
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(iv) The XGobi program | ]] for displaying multi-dimensional data was used to help visualize 
the cells and their neighbors. 

To compute the volumes and second moments of the cells we decomposed the cells into simplices 
and used the formulae in [|7J and [11, Chap. 21]. 



2. Notation 

We assume the n-dimensional lattice A is the intersection of n-dimensional lattices Ai, . . . , A r . 
The intersections of the Voronoi cells for the Aj partition IR n into a tesselation or honeycomb. Let 
V\ , . . . , Vk be representatives for the different polytopes or cells that appear in the honeycomb. 

Let pi (i = 1, . . . ,k) be the probability that a randomly chosen point in M. n (uniformly dis- 
tributed over a very large ball, say) belongs to Vi, and let iVj = piV/Vi, where Vi is the volume 



of Vi and V = \J det A is the volume of a fundamental region or Voronoi cell for A. Let V be the 
particular Voronoi cell for A that contains the origin. Then the honeycomb is periodic with "tile" 
equal to the part lying in V, and there are N\ cells of type V\ per copy of V, N2 cells of type V2, 
etc. Also 

V = N 1 V 1 + ■ ■ ■ + N k V k . (4) 

We assume that a point that falls into a cell of type Vi is quantized as the centroid c, of that 
cell, in which case the mean squared error is 

Ui = i \\x — Ci\\ 2 dx . 
JVi 

The mean squared error for the full quantizer is then 

i=i i=i 

We now derive the expression that we will use as a measure of the normalized mean squared 
error of this quantizer. Suppose we are quantizing a random variable X S R n , with differential 
entropy per dimension h(X), and whose support contains a large number of points of A. 

A general theorem of Zador about vector quantizers ( [j34|| , p5| , fli~8| , Eq. (20)]) implies that in 
this "high-rate" case the average mean squared error per dimension U/(nV) can be approximated 
by 

R_^ G2 2{h { X)-R) j (6) 

where R bits/symbol is the quantizing rate and G depends on the positions of the quantizing points 
but is independent of X. 
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Since G does not depend on the distribution of X, we may choose any convenient distribution 
in order to calculate G, and we assume that X is uniformly distributed over a large region of W 1 , 
or, equivalently, that X is uniformly distributed over V. Then 

h(X) = - log 2 V, so 2 2h ^ = V z/n . ( 7 ) 



n 



To calculate R, observe that we need h(pi, . . . ,pk) = — Yli=i Pi ^°S2Pi bits to specify the type 
of cell to which the quantized point belongs, and a further X^=i Pi 1°S2 -^j = ^2i=iPi^°S2(pi^/^i) 
bits to specify the particular one of the iVj cells of that type. This requires a total of log 2 V — 

Y^i=x Pi 1°S2 V% bits, and then R is this quantity divided by n, so that 

k 

2~ 2R = V~ 2 / n Yl Vf n/n . (8) 
k=i 

From (|)-(|) we obtain 

G= V ' a , (9) 



the normalized mean squared error per dimension, which we take as our figure of merit for a 
quantizer. The numerator of (^) is equal to U/V (see (|5|)). Note that when there is only one kind 
of cell (||) reduces to the familiar formula 

G = TT~2" ( 10 ) 

nV 1+ n 

for a lattice quantizer (@, [Iljl , [16]). 

Incidentally, a different expression from (^) for the figure of merit was used in a recent paper 



of Kashyap and Neuhoff |2Cj. By defining the rate of the quantizer in a different way they end 



up with a different expression in the denominator. However, we believe our formula gives a fairer 



comparison. This point is discussed in more detail in [27|. 



Further terminology. Let A be a lattice in M. n . The dual lattice will be denoted by A*. The 
norm of a vector x € M. n is its squared length (x, x). A similarity a is a linear map from W 1 to M. n 
such that there is a real number with (ax, ay) = N(x,y) for x,y £ M n . If A and M are similar 
lattices we write A = M. A lattice A is said to be N -modular if A = A* under a similarity that 
multiplies norms by N (cf. [ 5[). For example, the root lattices Z n and E% are 1-modular, A2 is 
3-modular, and D4 is 2- modular. We write ( Ai, A 2 , . . . ) for the lattice generated by the union of 

the lattices Ai, A 2 , Two lattices or polytopes are congruent if one can be mapped to the other 

by an element of SO(n). 



3. Writing a lattice as an intersection 



Let A be a lattice in R n . We wish to write A as an intersection 

A = Ai n A 2 n ••• n A r (11) 

where r is small and the Aj are pairwise congruent and as "simple" as possible. Ideally we would 
like each A.; to be a direct sum of congruent copies of a fixed low-dimensional lattice K such as Z, 
A2 or D4, but this is not always possible. In the example shown in Fig. |], for instance, the Aj are 
rectangular rather than square lattices. If this is not possible, we ask that the Aj be decomposable 
into a direct sum of as many congruent low-dimensional sublattices as possible. 

If we were going to investigate the honeycombs associated with higher-dimensional intersections 
such as those for Eq or Eg, we would impose an additional formal requirement that the Aj form 
an orbit under some subgroup of the automorphism group Aut(A), in order to guarantee that the 
honeycomb be symmetric. However, for the low-dimensional examples studied in Sections 5-7, we 
were able to achieve this symmetry by using the natural decompositions, without introducing the 
machinery of group theory. 

This section describes some general methods for finding intersections. We will make use of the 
following version of de Morgan's law: 

Lemma 1. If A±, . . . , A r are lattices in W 1 then 

Ain---nA r = ( A^,...,A* }* . 



Proof. See for example [24, Section 82F] 



Method 1: Partitioning the minimal vectors for A*. For the first method, suppose the 
dual lattice A* is generated by its minimal vectors (i.e. the vectors of minimal nonzero norm), and 
let K be one of Z, A2 or D4. Let K have dimension k and be A^-modular. Then we attempt to 
partition the minimal vectors of A* into a number of copies of the minimal vectors of a suitably 
rescaled version of K n / K . 

If this can be done, let A* be the lattice generated by the ith part of the partition. The 
hypotheses guarantee that A* = Aj = K n / K , and, by the Lemma, 

A = Ai n • • • n A r , 
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where r = number of minimal vectors of A* divided by number of minimal vectors of K. 

Whether this partitioning is possible is an interesting question in its own right. For example, 
can the 240 minimal vectors of Eg be partitioned into 15 copies of the minimal vectors of (a scaled 
version of) Z 8 , i.e. into 15 coordinate frames^] or into 10 copies of the minimal vectors of A^l 
Partial answers are given below. 

We now give a number of examples, beginning with the case when K = Z. 



D4. The lattice may be taken to have generator matrix 



2 
2 
2 
1111 



(12) 



and then the 24 minimal vectors consist of eight of the form (±2, 0, 0, 0) and 16 of the form 
(±1, ±1, ±1, ±1). D4 is 2-modular and is generated by its minimal vectors. The minimal vectors 
may be partitioned into three coordinate frames, consisting of ±1 times the rows of each of the 
matrices 



2 

2 

2 

2 



+1 +1 +1 +1 

+1 -1 +1 -1 

+1 -1 -1 +1 

+1 +1 -1 -1 



-1 +1 +1 +1 

-1 -1 +1 -1 

-1 -1 -1 +1 

-1 +1 -1 -1 



(13) 



After applying the Lemma and rescaling, we conclude that if Ai, A2, A3 (= Z 4 ) have the generator 

Ai n A 2 n A 3 = A , 



matrices given in (13), then 



where A has generator matrix 



4 

2 2 

2 2 

2 2 



(14) 



(15) 



and is another version of D4 on the scale at which its minimal norm is 8. Equation (14) may also 
be verified directly, without appealing to the Lemma. The group generated by the second matrix 
in (|l~3l) and diag{ — 1, +1, +1, +1} is a symmetric group S3 permuting the A$; it is also a subgroup 
of Aut(A). The honeycomb for this example is studied in Section 7. 



Barnes- Wall lattices. The preceding example can be generalized using orthogonal spreads. Let 
BW n (n = 2 m 



1,2,...) denote the n-dimensional Barnes- Wall lattice ( |llf| , |21[| , J 2j). In 

±v„ with (vi,Vi) = a constant, (vi,Vj) = if i 7^ j. 



§ A coordinate frame in R n is a set of 2n vectors ±«i, 
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particular, BW 2 = Z 2 , BW4 = At, BWg = Eg. For n / 8, BW n is 2-modular, while as already 
mentioned BWg is 1-modular. 



It is known that the minimal vectors of BW n may be partitioned into f^^I^ 1 (2 J — 1) coordinate 
frames, which are transitively permuted by symmetries of Aut(BW n ). This is a consequence of the 
existence of an orthogonal spread in the orthogonal vector space £l + (2m, 2) of maximal Witt index 

(§, §, 0)- 

It follows that BW n can be written as the intersection of nj=i 1 (^"' — -0 c °pi es °f 111 
particular, Eg is the intersection of 15 copies of Z 8 . An explicit method for constructing such an 
intersection for Eg is given below. 

Intersections of smaller numbers of lattices are possible, although they are less symmetric and 
therefore less satisfactory. For example in ( |i~4| ) it is also true that Ai D A2 = A = D4. Similarly, Eg 
is (up to a similarity) the intersection of 2Z 8 and the lattice (similar to 2Z 8 ) with generator matrix 



1 


-1 



1 



1 



1 




1 

-1 




-1 



1 




-1 



-1 
-1 







1 



1 
1 



-1 







1 



-1 



1 



-1 



But this representation of Eg is in no way canonical, and the resulting honeycomb is not interesting. 



Eisenstein and Hurwitzian lattices. Smaller intersections which are canonical can be obtained 
if we change K from Z to A 2 or D4. For example, the minimal vectors of Eg can be partitioned 
into 10 copies of the minimal vectors of A\. As in 11], let £ = {a + buj : a, b £ Z}, uj = e 2m ^ 3 , 



denote the ring of Eisenstein integers. The six units in £ are ±1, ±uj, ±cj. When regarded as a 
two-dimensional real lattice £ is similar to Ai- As an £ -module, Eg has generator matrix 



6 

e 

1 1 

1 







1 

-1 1 



where 6 



to — U) . 



(16) 



Inner products are computed using the Hermitian inner product (u, v ) = ^ U{Vi. See [U, Chapters 
2 and 7] for further details. The minimal vectors consist of 24 of the form (u9, 0, 0, 0), where u is 
a unit in £, and 8 x 3 3 = 216 which are congruent mod 9 to one of the eight nonzero codewords of 
the tetracode |ll], Chap. 3]. 
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A partition of these 240 vectors into 10 copies of the minimal vectors of A\ was found by 
graph coloring. A graph was constructed with the 40 projectively distinct vectors as nodes and 
with edges corresponding to pairs of non-orthogonal vectors. A coloring with 10 colors was then 
found with the help of a program supplied by David Johnson. The ten copies of £ 4 = A\ are 
shown in Table |2[ Only one from each complex conjugate pair is shown. By applying the Lemma 
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— LO 
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— LO 





oo 
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Table 2: Decomposition of minimal vectors of Eg into ten copies of £ 4 = A\. Each row generates 
a copy of £ = A|. The complex conjugates of the last four rows have been omitted. 

we obtain a representation of Eg as an intersection of 10 copies of A\. In fact (since £ is itself 
3-modular) we may omit the final step of taking the duals of the lattices in Table ^. Let Ai, . . . , Aio 
be the ten versions of A\ generated by the rows of Table || and their complex conjugates. Then 
their intersection is easily seen to be the version of Eg with generator matrix times (|l6|). This 
decomposition is probably not unique, and it would be nice to know which version has the largest 
symmetry group. 

Again just two lattices suffice: E$ is also the intersection of the first two lattices in Table ^. 
We may also write Eg as the intersection of five copies of D\. For this we regard Eg as a 



2-dimensional module over the ring 
of D\ have generator matrices 



Z?4 of Hurwitzian quaternions [11, p. 55]. The five copies 



1 + t 
1 + i 



k 
-k 



(17) 



Ei may be written as the intersection of four copies of A\ with generator matrices 
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1 


00 


LO 




1 


u) 


LO 




1 


LO 


1 _ 



(18) 



This was found by partitioning the 72 minimal vectors of Eq (by hand) into four copies of the 
minimal vectors of A| and using the Lemma. 



Method 2: Congruence bases and norm-doubling maps. The second method is based on 
the observation that several well-known lattices A have the property that for some prime ir, the 
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vectors in some of the classes of A/nA can be partitioned into coordinate frames. For example, 
Conway's proof of the uniqueness of the Leech lattice A24 Jll], Chap. 12] considers the classes of 
A24/2A24. A consequence of the numerical identity 

n n 4 n 6 n 8 1 196560 16773120 398034000 

— + — + — + — = — I 1 1 

1 2 2 48 1 2 2 48 

= 16777216 = 2 24 , (19) 

where rij is the number of vectors in A24 of norm j, is that, for the classes of A24/2A24 in which 
the minimal norm is 8, the minimal vectors in the class form a coordinate frame or congruence 
base. A similar property holds for the _D 4 , E$, K\2 and other lattices (cf. This gives a 

representation of A24 as an intersection of 398034000/48 = 8292375 copies of Z 24 . However, the 
following argument, due to J. H. Conway (personal communication), shows that if the lattice has 
a suitable norm-doubling map (cf. jll], p. 239], Q) then we can also partition the minimal vectors 
into coordinate frames and obtain a smaller intersection. 

Suppose a lattice A C W 1 has the structure of a free module over a ring J with inner product 
( , ) (cf. pi] , Chap. 2]). In the present application J will be either Z or £ . Let a = n/dimjjj. 
Consider the classes of A/2A. Suppose there is an integer m with the property that each class either 
contains no vectors of norm 2m, or else all the vectors of norm 2m in the class can be partitioned 
into sets of 2a vectors ±«i, ±t>2, • • • , where Vj) = Qiii^=j. 

Suppose in addition there is a norm-doubling map T, a similarity from A into A such that 
(Tu,Tu) = 2(u, u) for u £ A, with the extra property that 2A C TA. Then we may conclude that 
the vectors of norm m in A may also be partitioned into sets of 2a mutually orthogonal vectors. 

To see this, let u € A have norm m. Then v = Tu has norm 2m, and by the hypotheses is part 
of a coordinate frame ±t?i, . . . =t v a , where Vi = v + 2u>j, say, with w\ = 0. We can write 2wi = Tw[ 
for some w[, so Vi = T{u + w^). Since T is a similarity, the set ±(u + w'j) is a coordinate frame 
containing u. 

Examples 

(i) A = D4 or J = Z, m = 2, T = direct sum of respectively 2 or 4 copies of , 
with T 2 = 21. The classes of -D4/2D4 and E 8 /2Eg and the associated congruence bases are given 
in [jO], Chap. 6]. The analogue of ( fl9|) for Eg reads 

240 2160 o 

H + = 2 8 . 

2 16 
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We obtain decompositions of the minimal vectors of -D4 into three coordinate frames, as already 
seen in (|l3|), and of the minimal vectors of E% into 15 coordinate frames as also discussed above. 
To get an explicit decomposition in the latter case, note that a coordinate frame of norm 4 vectors 
has the form Tu + 2wi = T(u + Twi). So the coordinate frame of norm 2 vectors consists of the 
vectors of minimal norm in the translate u + TA. For Eg these consist of seven sets of the form 
shown on the left in (GO) and eight of the form shown on the right: 
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(20) 



Then Eg is also the intersection of the 15 copies of Z 8 having these generator matrices. 

(ii) A = Leech lattice A24, J = Z, m = 4, T = (I + i), where i G ^4ui(A24) is given in [11 
Fig. 6.7], and satisfies i 2 = —I, with (/ + — i) = T(I — i) = 21. The coordinates frames 
of norm 4 vectors consist of the 48 vectors of minimal norm in the translates u + (I — i)A where 
(u, u) = 4. We obtain a decomposition of the minimal vectors of A24 into 4095 coordinate frames, 
and a representation of A24 as the intersection of 4095 copies of Z 24 . 

We do not know if it is possible to write the Leech lattice as the intersection of 2730 copies of 
A^ 2 . Since 196560/96 is not an integer, there is no analogous decomposition as an intersection of 
copies of D\. 

(iii) A = Coxeter-Todd lattice K12, a 6-dimensional ^-module, J = £, m = 6, T = the map 
given in |], Eq. (43)], with T 2 + T + 2 = 0. The classes of K\2j2K\2 are given in ||, and the 
analogue of (|l9l) reads 



756 4032 20412 

H 1 1 

2 1 12 



,12 



The coordinate frames of norm 6 vectors consist of the 12 minimal vectors in the translates 
u + (T + I)Ki2, (u,u) = 6. By combining these coordinate frames in sets of three, by taking 
the union of the sets a{±vi, . . . ,±^12} with a = l,u> and Q, we obtain a decomposition of the 
minimal vectors of K\2 into 21 copies of the minimal vectors of A®, and, via the Lemma, a repre- 
sentation of K12 as the intersection of 21 copies of A\. An explicit decomposition, not shown here, 
was found by the graph coloring method mentioned earlier. 
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Method 3: First principles. If the above methods fail, as they do for the bcc and fee lattices, 
we can always fall back on a direct attack from first principles. The following method handles 
the hexagonal, bcc and fee lattices in a unified manner. We list the vectors of small norms in the 
lattice, and look for a partition of some subset of these vectors which produces a small number of 
congruent, decomposable lattices whose intersection is similar to the original lattice. 

For the hexagonal lattice, which we take to be generated by (0, 1) and (— ^), there are six 
vectors of norm 1, namely (0,±1), (±-^,±^), and six of norm 3, namely (±\/3, 0), (±^,±|). 



Then (||) is obtained by partitioning these 12 vectors into three sets of size 4. 

For th 
following: 



For the bcc lattice generated by (1, 0, 0), (0, 1,0), (^, ^, ^), the vectors of small norms are the 



shape 


ill 

2 2 2 

1 



norm 


3 
1 
1 

2 



number 

1 

8 
6 

12 



1 1 

We take the 18 vectors of norms 1 and 2 and partition them into three sets of size 6. The resulting 
lattices have generator matrices 



1 
1 
1 



1 

1 
1 




1 
1 



1 

1 
-1 



(21) 



and their intersection has generator matrix 



2 

2 

1 1 1 



(22) 



which is indeed another version of the bcc lattice. This decomposition of the bcc lattice is the 
simplest we have found, and will be discussed further in Section 5. 

For the fee lattice generated by (1, 1, 0), (1, —1,0), (0, 1, —1), the vectors of small norms are: 



shape 


norm 


number 








1 


1 1 


2 


12 


2 


4 


6 


2 1 1 


C 


24 


2 2 


8 


12 


3 1 


10 


24 


2 2 2 


12 


8 



The simplest intersection we have found is formed by taking the 32 vectors of norms 6 and 12 and 
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partitioning them into four sets of size 8. The resulting lattices have generator matrices 



2 1 
1 2 
-2 2 



2 1 
1 2 

-2 2 



2 
1 

-2 



(23) 



and their intersection is the fee lattice with generator matrix (3, 3, 0), (3, —3, 0), (0, 3, —3). We will 
return to this example in Section 6. 



Method 4: Using intersections of codes. A fourth method, which however has not yet led to 

any interesting examples, is to reduce the problem to the analogous question for codes. Let A(C) 
denote the lattice obtained by applying Construction A to a binary linear code C ( fill , Chap. 5]). 
If Ci, . . . , C r are codes of length n whose intersection is a code C, then 

A(c) = A(Ci) n • • • n A(c r ) . 



This can be generalized to nonbinary codes [11, Chaps. 7, 8], in particular to the case where the 
Q are nonbinary codes whose intersection is binary. 

4. The hexagonal lattice as an intersection of three lattices 

We now begin our study of the honeycombs formed by some of the intersections described in 
Section 3. Once the cells have been found it is generally straightforward to compute their volumes 
and second moments using the techniques presented in ]7j or [jll], Chap. 21]. 

The hexagonal lattice is the intersection of the three rectangular lattices given by (^), as in 
Fig. U The honeycomb is shown in Fig. |2[ We now compute the mean squared error for this 
quantizer. This analysis differs from that given in Q, which used a different set of representation 
points for the quantizer. 

There are four types of cells. The origin is contained in a hexagon V\ (horizontally shaded in 
Fig. |2|) of edge length l/\/3, area V\ = \/3/2 and second moment U\ = 5\/3/72. The second type 
of cell is a small equilateral triangle V2 (cross- hatched), with V 2 = y/3/12, U 2 = v / 3/432. V 3 is an 
isosceles triangle (diagonally shaded), with V3 = v3/12, U3 = 5v3/1296. The fourth type, V4, is 
a larger equilateral triangle (vertically shaded), with V4 = v3/4, U4 = y/3/48. 

The incidences between the different types of cells are shown in Fig. [3|. Here (and in similar 
diagrams in later sections) a circle containing i refers to a cell of type Vi, and an edge 
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indicates that V% and Vj share a common face of the maximal dimension (here 1), and the edge 
labels indicate that each V% is adjacent to s cells of type Vj and each Vj to t cells of type V%. 




Figure 3: Incidences between cells in honeycomb for hexagonal lattice. 

The Voronoi cell V for the intersection lattice is enclosed by the broken lines in Fig. ^. In the 
notation of Section 2, V contains N\ = 1 copy of Vx, N2 = 6 copies of V2, N3 = 6 copies of V3 and 
= 6 x i = 2 copies of V4, and the volume equation (||) reads 

^f + f + f + f- <*> 

Thus the probabilities pi , . . . , P4 of a randomly chosen point in the plane belonging to a cell of each 
type are all equal to 1/4. From the mean squared error for this quantizer is G = 2 7 / 4 /27 = 
0.1246 . . . This value (and corresponding values of G found in the next three Sections) is considerably 
worse than the value 0.080188 . . . for the hexagonal lattice itself. 

5. The bcc lattice as an intersection of three lattices 

The bcc lattice is the intersection of the three "rectangular" lattices Ai, A2, A3 defined in (^Tj). 
Each Aj is congruent toZx x y/21i, and has as Voronoi cell a brick with square cross-section. 
There is an obvious symmetric group £3 that permutes the Aj. 

Again the honeycomb contains four types of cells. The origin is contained in the intersection of 
the Voronoi cells at for the three Aj. This is a cube, Vi, with vertices (±1/2, ±1/2, ±1/2), volume 
V\ = 1, and second moment U\ = 1/4. Across each square face of V\ is a square pyramid V2, such 
as that with base (1/2, ±1/2, ±1/2), apex (1,0,0), V2 = 1/6, U2 = 11/600. Across each triangular 
face of V2 is a tetrahedron V3, such as that with vertices (1/2, ±1/2, 1/2), (1,0,0), (1,0,1/2), 
V3 = 1/24, U3 = 1/512. Finally, across the other three faces of V3 we reach a fourth type of cell, 
V4, a quarter-octahedron, which occurs in two orientations, one having vertices such as 

^,0,1), (1,0, l), (l,0, i), (1,0,1), Q,44), (25) 

the other having vertices such as 

i- 1) ■ M)- (i.|> )' ^^^(il'D^IM) ■ < 26 > 
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These may be described as quarters of squat octahedra. E.g., is a quarter of the octahedron 
with vertices (±1,0, ±1), (1/2, ±1/2, 1/2). For T 4 we have V A = 1/12, U 4 = 1/192. 

No further types of cell appear: every face of either version of Va leads to a Vs- The incidence 
diagram is shown in Fig. [|. 




Figure 4: Incidences between cells for honeycomb for bcc lattice. 

The Voronoi cell V for the intersection lattice is a truncated octahedron with 24 vertices 
(0, ±1/2, ±1). This contains N± = 1 copy of Vi, N2 = 6 copies of V2 and N3 = 24 copies of 
P3. The cells of type V4 partially overlap V. There are 12 of type (p5|), intersecting V in a tetra- 
hedron such as that with vertices (1/2, ±1/2, 1/2), (1/2,0,1), (1,0, 1/2), with volume 1/24. There 
are also 24 of type (|26|), intersecting V in a tetrahedron such as (1/2,1/2,1/2), (1,0,0), (1,0,1/2), 
(1, 1/2,0), with volume 1/48. The volume equation (||) reads 

1111 

4 = lxl + 6x- + 24x h 12 x h 24 x — 

6 24 24 48 

= 1 + 1 + 1 + 1 , (27) 

so again the probabilities pi of a randomly chosen point belonging to a cell of given type are all 
equal to 1/4. The mean squared error is 

G. 151^=0.1366... 
9600 

6. The fee lattice as an intersection of four lattices 



The fee lattice is the intersection of the four lattices Ai , . . . , A4 defined in ( f23[) . Each of these 
has Gram matrix equivalent to 



6 
-3 




-3 
6 
12 



and is a direct sum y/SA2 © \/l2Z, with Voronoi cell a hexagonal prism. The Aj look more 
symmetrical if they are written in the coordinates used to describe the root lattice A3 (= D3), 
that is, using four coordinates that add to 0. Then Ai has generator matrix 






3 
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Figure 5: Incidences among cells of fee honeycomb, 
and the others are given by cyclic shifts of these columns. This shows that there is a symmetric 



group S4 permuting the Aj. However, the three-dimensional coordinates given in (23) are more 
convenient for computations. 

This honeycomb is the most complicated we have analyzed and we shall give only a brief 
description. There are twelve types of cells, V\, ■ ■ ■ ,Vi2, whose parameters are summarized in 
Table [3| and whose incidences are shown in Fig. || Figures || and shows cross-sections through 
the honeycomb along the planes z = and z = 0.35. 

The following is a brief description of the cells, including coordinates for one cell of each type. 

V\. Obtained from cube by pushing in corners and pulling out edge midpoints. Vertices: all 
cyclic shifts and sign changes of (3/2,0,0), (1,1,0), (3/4,3/4,3/4). 
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i 


V 


e 


J 


A r 

Ni 


T T 

Vi 


Pi 


T T 

Ui 


1 


26 


48 


24 


1 


9 


1/6 


1449/160 


2 


5 


8 


5 


24 


1/8 


1/18 


41/3200 


3 


6 


12 


8 


24 


1/4 


1/9 


9/320 


4 


5 


8 


5 


24 


1/8 


1/18 


41/3200 


5 


5 


9 


6 


24 


1/10 


2/45 


427/50000 


6 


4 


6 


4 


48 


1/40 


1/45 


443/320000 


7 


10 


18 


10 


8 


27/40 


1/10 


41013/160000 


8 


6 


12 


8 


6 


1 


1/9 


47/180 


9 


5 


9 


6 


8 


3/8 


1/18 


189/3200 


10 


4 


6 


4 


24 


1/40 


1/90 


1897/1280000 


11 


5 


8 


5 


24 


3/40 


1/30 


2319/400000 


12 


22 


36 


16 


2 


63/10 


7/30 


213597/40000 



Table 3: The twelve types of cells in the fee honeycomb, showing numbers of vertices, edges, 
faces (v,e,f), the number per fee cell (Ni), and their volumes, probabilities and second 
moments(Vi, pi, Ui). 

P 2 . Pyramid with kite-shaped base. Base: (3/2,0,0), (1,1,0), (1,0,1), (3/4,3/4,3/4), apex: 
(3/2,1/2,1/2). 

P 3 . Irregular octahedron with vertices (3/2,0,0), (3/2,3/2,0), (1,1,0), (2,1,0), (3/2, 1/2, ±1/2). 
V±. Congruent to P 2 . Vertices: (3/2,0,0), (3/4,3/4,3/4), (3/2,1/2,1/2), (1/2,3/2,1/2), 
(1,1,0). 

P 5 . Irregular polyhedron with six faces. Vertices: (3/2,0,0), (2,1,0), (2,0,1), (9/5,3/5,3/5), 
(3/2,1/2,1/2). 

P 6 . Tetrahedron: (2,1,0), (3/2,3/2,0), (3/2,1/2,1/2), (9/5,3/5,3/5). 

V'j. "Flying saucer": hexagonal base (cyclic shifts of (3/2,3/2,0) and (9/5,3/5,3/5)) with four 
vertices above it (cyclic shifts of (3/2,1/2,1/2) and (3/4,3/4,3/4) at apex). 

V 8 . Irregular octahedron with vertices (2, ±1,0), (2,0,±1), (3/2,0,0), (3,0,0). 

Pg. Regular tetrahedron (vertices (2,1,0), (2,0,1), (3,0,0), (3,1,1)) with triangular cap (apex 
(9/4,3/4,3/4)) on one face. 

Pio. Irregular tetrahedron with vertices (2,1,0), (2,0,1), (9/4,3/4,3/4), (9/5,3/5,3/5). 

P X \. Another pyramid with kite-shaped base. Base: (3/2,3/2,0), (9/4,3/4,3/4), (9/5,3/5,3/5), 
(12/5,6/5,3/5), apex: (2,1,0). 

P\2- See Fig. |8[ Has 26 vertices, four hexagonal and 12 kite-shaped faces. Vertices are all permu- 
tations of (3/2,3/2,0), (3/2,3/2,3), (9/4,3/4,3/4), (9/4,9/4,9/4), (9/5,3/5,3/5), (12/5,9/5,9/5), 
(12/5,6/5,3/5). 
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Figure 6: Cross-section of fee honeycomb along plane z = (—3 < x < 9, —6 < y < 6), with origin 
at center of octagon on left. Only three cells are visible: V\ (octagon), V2 (small kite), "Pg (large 
dart). 

If we decompose M 3 into Voronoi cells for the intersection lattice A (24), just three of the twelve 
types of cells are cut by the boundary walls. Cells of type Vg are cut into three equal pieces, cells of 
type Vn are cut in half, and cells of type V12 are cut into four equal ice-cream cone shaped pieces. 
The base of each cone is at the center of V12 , and the top contains one of the hexagons and parts 
of the neighboring faces. 

The mean squared error is 

= 12269777 28/135 8/27 38/45 = m2 
816480000 

7. The D4 lattice as an intersection of three cubic lattices 

This is the nicest example we have found. The three cubic lattices Ai, A 2 , A3 are defined in 
(|l2|) and their intersection A = D4 in Jl3|). There are just four types of cells, V\, V2, ^3, 'Pi, 
whose properties are summarized in Table [| and whose intersections are shown in Fig. ^. We use 
coordinates (a,b,c,d) for points in M 4 . 

V\ is the 4-dimensional regular polytope known as a 24-cell ( pl[ , Jl2[). The Voronoi cells for 
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Figure 7: Cross-section through fee honeycomb along plane z = 0.35. Cross-sections of all 12 types 
of cells can be seen. 

Ai, A2, A3 at the origin are all cubes, whose intersection is bounded by the hyperplanes 
\a\ < 1, \b\ < 1, |c| < 1, \d\ < 1, \a\ + |b| + |c| + \d\ < 1, 

which is the 24-cell with vertices of the form (±1, ±1,0,0). 

Across each of the 24 octahedral faces of V\ we reach an octahedral-based pyramid V2, such as 
that with vertices (1, ±1,0,0), (1,0, ±1,0), (1,0,0,±1) and (2,0,0,0) (the apex). 

There are eight other faces of V2, regular tetrahedra; these lead to copies of V3, which is an 
irregular simplex such as that with vertices (1,1,0,0), (1,0,1,0), (1,0,0,1), (2,0,0,0), (1,1,1,1). 

i v f Ni Vi pi Ui 

1 24 24 1 8 1/4 104/15 

2 7 9 24 1/3 1/4 8/105 

3 5 5 96 1/12 1/4 11/900 

4 6 9 32 1/4 1/4 1/20 

Table 4: Cells in D4 honeycomb, showing numbers of vertices and 3-dimensional faces (v,f), the 
number per D4 cell (Ni) and their volumes, probabilities and second moments (Vi,pi,Ui). 



21 



Figure 8: 26- vertex cell V\i- 




Figure 9: Incidences among cells of -D4 honeycomb. 



Finally, three of the five faces of each V3 lead to cells of the fourth type, V&. This can best be 
described as the product of two skew equilateral triangles of different sizes (just as a tetrahedron 
in three dimensions is the product of two skew line segments). Take an equilateral triangle with 
vertices p = (2,1,1,0), q = (1,1,0,0), r = (1,0,1,0) and another with vertices P = (2,0,0,0), 
Q = (1, 1, 1, 1), R = (1, 1, 1, —1). Then is their convex hull. There are nine tetrahedral faces 
given by the convex hull of an edge of the first triangle and an edge of the second triangle. 

If we decompose M 4 into Voronoi cells for the intersection lattice, only cells of type V\ are cut by 
the boundary walls. Each V4 is divided into three equal pieces, a typical piece being an "ice-cream 
cone" whose center is at the center of V4 and whose three-dimensional face is the convex hull of 
the second triangle (P, Q, R) and any edge of the first triangle. 
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The volume equation (|]) then reads 



32 = lx8 + 24 x- + 96 x — 




Again a random point is equally likely to fall into a cell of any of the four types. The mean squared 



8. Conclusions and comments 

In each of the honeycombs of Sections 4, 5 and 7 just four types of cells occurred. This is easily 
explained in the case of the D4 honeycomb: there are three equivalent lattices Ai, A2, A3, and the 
associated quantizers are essentially making binary decisions about the location of a point with 
respect to the intersection lattice. So the space is divided up into regions that can be labeled 000, 
001, 011 and 111. 

This argument doesn't quite apply to the A2 or bcc honeycombs, since there the individual 
lattices themselves are not fully symmetric (rectangular rather square in the A2 case, for instance). 
So it is fortuitous that only four cells occur. 

In contrast, the fee honeycomb of Section 6 shows that the number of cells can increase rapidly 
in less fortunate cases with more component lattices. It would be interesting to see how complicated 
are the Eg and Eg honeycombs mentioned in Section 3. 

As for practical applications, it is possible that better quantizers could be obtained by amalga- 
mating less symmetrical cells. For example, in Fig. g the diagonally and vertically shaded triangles 
could be amalgamated to give a honeycomb made up of regular hexagons and equilateral triangles 
with the same edge length as the hexagons. The new honeycomb will have larger absolute error 
but a smaller normalized error G. We did not investigate this possibility for the higher-dimensional 
examples. 

Another topic for future research is the design of quantizers by taking a simple initial lattice 
Ai and combining it with several other lattices which are both rotations and translations of Ai. 
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